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ABSTRACT 

Aims. We investigate observations of the X-ray binary Cygnus X-1 with unusually high hardness and low flux. In particular, we study 
the characteristic frequencies seen in the PDS and the hardness-flux correlation within and between these observations. 
Methods. We analyse observations of Cyg X-1 during periods when the source reaches its highest hardness levels (j£ 1 for the 9- 
20 keV over 2-4 keV RXTE/PCA count ratios, corresponding to T ^ 1.6). Using the relativistic precession model to interpret the PDS, 
we estimate a value for the inner radius of the accretion disc. We also study the hardness-flux correlation. 

Results. In the selected observations, the characteristic frequencies seen in the power spectrum are shifted to the lowest end of their 
frequency range. Within a single observation, the hardness-flux correlation is very weak, contrary to the negative correlation normally 
observed in the hard state. We suggest that this could be interpreted as the inner disc boundary being at large radii 50 R g ), thereby 
requiring more time to adjust to a changing accretion rate than allowed by a single RXTE observation, and compare our findings to 
estimates of the viscous time scale responsible for small scale variability in the system. 
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tional Comptonized component, is also required to accurately fit 
the spectrum. 

Many sources also show a clear difference between the hard 
and soft states in the power density spectrum (PDS). In Cyg X-1, 
the PDS in the soft state is generally well described by a sin- 
gle p ower law, with a cut-off at higher frequencies dCui et al.l 
1 19971) . while the PDS of the hard state shows a more com- 
plex structur e, with several breaks or quasi-periodic os cilla- 
tions (QPOs; iBelloni & Hasingeril990tlNowaketal.ll999l) . The 
hard state PDS can be gen erally well-fit with broad Lorentzian 
components (Nowak 2000). The number of components varies 
between studi es, and depends on the frequency window and 
models used. Pottschmi dt et al. I (120031) conducted a large study 
using up to four Lorentzians in the 0.002-128 Hz range. In 
Axelsson, Bor gonovo & Lars son (2005]) it was found that two 
Lorentzians are generally sufficient when fitting the PDS in the 
0.01-25 Hz range; by adding a cut-off power-law component to 
the model, the complete evolution of the PDS from hard to soft 
state could be modelled. 

In this paper, we analyse observations of Cyg X-1 in obser- 
vations with unusually high hardness and low flux, the "hard 
edge" of the hard state. In particular, we study the characteris- 
tic frequencies seen in the PDS and the hardness-flux correlation 
within and between these observations. 

2. Observations and data reduction 

Since i ts launch on Dece mber 30, 1995, the All-Sky Monitor 
(ASM, iLevine et all 1 19961) aboard the RXTE satellite has pro- 
vided nearly continual coverage of the X-ray sky in the 2-12 keV 
range. A given source will be observed several times each day, 
and the data is available either in dwell-by-dwell format, or one- 



1. Introduction 

Cygnus X-1 is one of the most studied X-ray sources, and is 
often quoted as the prototype black hole binary system. The 
source exhibits two main spectral states, commonly referred to 
as hard and soft, with a brief intermediate state du r ing transi- 
tions (e.g.. iTanaka & Lewinlll995l; ICui et al.|[l997t lEsin et alj 
119971; [Zdziarski & Gierlinskil l 2004 . Several models have been 
proposed to explain the observed states and transitions. The 
two main components of such models are usually a geomet- 
rically thin, optically thick accretion disk and a geometrically 
thick, optically thin hot inner flow or corona. The models vary 
in the geometry and properties of mainly the hot flow/c orona 
dHaardt & Maraschill 19931; iBeloborodovl 19991; ICoppill 19991) . 

Transitions between the two states are believed to be a 
response to a change in accretion rate. In one widely ac- 
cepted model (see e.g.. iGierlihski et al.ll 19971; lEsin et all [l997; 
Gierlihski et al. 1999llZdziarski et al.ll2002l7 these changes result 



in a reconfiguration of the accretion flow. In the soft state (high 
accretion rate) the disc extends almost in to the last stable orbit, 
whereas in the hard state (low accretion rate) it is truncated and 
replaced by a hot flow at inner radii. Throughout the hard state 
and into the transition, the hardness of the spectrum and chang- 
ing temperature of the observed disc blackbody component is 
believed to track the variable radius of the inner disc. 

The ener gy spectrum o f Cyg X-1 in its hard state was de- 
scribed by iGierlinski et al.l ((T997). It is dominated by a com- 
ponent arising from thermal Comptonization in a plasma with 
electron temperature of ~ 100 keV and optical depth t ~ 1; only 
a weak blackbody component is visible, corresponding to an in- 
ner disc temperature of ~ 200 eV in a geometry with a truncated 
disc. An additional soft component, often modelled as an addi- 
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Fig. 1. The one-day average ASM light curve for Cygnus X-l. The hard and soft states are clearly visible in the light curve, as well 
as a number of large flares. The vertical lines in the figure mark the hard state episodes with unusually low flux analysed in this 
paper. 



day averages. Figure Q] shows the one-day average lightcurve for 
Cygnus X-l. The periods when the source enters the soft state 
are clearly visible, as well as a number of flaring episodes. There 
are also a number of occasions when the ASM flux drops be- 
low the normal hard state level, and these are the observations 
studied in this paper. The data used are from pointed observa- 
tions made with the PCA instrument onboard the RXTE satellite 
dJahoda et al.l[l9 96). The times of the observations are indicated 
with vertical lines in Fig.Q] The typical length of an observation 
is ~ 0.5-1 hour. 

To create the PDS, lightcurves were extracted with 10 ms 
resolution in the 2-9 ke V band. The lightcurves were then sep- 
arated in intervals of 8192 bins, and a PDS was computed for 
each interval. Only intervals without gaps were retained, and the 
resulting PDS were then averaged to form the final power spec- 
trum. In calculating the hardness, PCA Standard2 data with 16 s 
resolution were used. As the data stem from many different ob- 
servations with varying setups, the PCUs in use vary between 
observations; all lightcurves have therefore been normalized to 
one PCU. The latest available bright source background model 
(dated August 6, 2006) was used. 

In the following analysis, we will concentrate on the obser- 
vations with lowest flux, which also have the highest hardness 
levels, in agreement with the negative ov erall flux-hardness cor - 
relation in the hard state (see for example Zdzia rski et al.l l2002). 
For easy reference, we will use the term "hard edge" to describe 
the episodes when the hardness ratio (HR) of Cyg X-l is roughly 
HR^; 1 with the HR calculated as the count ratio between the 9- 
20keV and 2—4 keV bands (corresponding to a spectral index 
r<1.6). 



3. Analysis and results 



Gilfa nov. Churazov & Revnivtsevl d 19991) found that the char- 
acteristic frequencies in the PDS of Cyg X-l correlated with 
the photon index . This was also se e n in t he sy stematic timing 
analys is of both iPottschmidt et al.1 d2003l) and lAxelsson et al.l 
(2005), where the peak frequencies of the Lorentzian compo- 
nents in Cyg X-l decreased with increasing hardness, indicat- 
ing that the frequencies are shifting systematically. Indeed, a 
third QPO-type component was seen to enter the upper region 
of the frequen c y win dow for the hardest observations in the 
lAxelsson et al.l d2005l) study, naturally explained as one of th e 
higher frequency components seen in IPottschmidt et al. (2003). 
Figure [2] shows a hardness-flux diagram of Cyg X-l in the hard 
state, with these observations represented by black squares. In 
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Fig. 2. Hardness-flux correlation in Cyg X-l during the hard 
state. Black squares indicate the locati on of the observations 
where a third Lorentzian was required by Axelsson et al. (2005) 
to fit the PDS. The hardness is calculated for the 9-20 keV over 
2—4 keV bands, and the flux is the sum of these two bands. 
Typical error bars are indicated in the lower left corner. 



all these observations, the two more common Lorentzian com- 
ponen ts are shifted to the low est part of their frequency range. 

In lAxelsson et al.l d2005l) . the two hard state Lorentzian com- 
ponents were identified with r elativistic precession frequencies 
dStella. Vietri & Morsinkl 19991) at the inner radius of a truncated 
disc. Shifts in frequency of the temporal components then corre- 
spond to changes in the truncation radius. Figure|2]clearly shows 
that a shift to low frequencies, characterised by the third com- 
ponent appearing in the frequency window, are connected with 
high hardness - all such observations occur in the hard edge of 
Cyg X- 1 . An example hard edge PDS, taken from an observation 
on MJD 52748 (points) and fit using Lorentzian components 
(lines), is seen in Figure[3] 

3. 1. Estimating the inner disc radius 

In the model used by Gierlinski et all dl997l) (eqpair, [Coppi 
1999) for the energy spectrum of Cyg X-l, a high hardness ratio 
reflects a high ratio of luminosity in the hot plasma electrons to 
that of the irradiating soft disc photons and a slightly lower disc 
temperature, suggesting a larger inner radius. The absolute val- 
ues of disc temperature and inner radius from spectral modelling 
of data from the RXTE are, however, not well constrained when 
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Fig. 3. Power density spectrum of Cyg X-l in the hard edge. 
Three Lorentzian components are necessary to fit the data. 



the peak of the disc blackbody, as here, lies well below the lower 
boundary of the instrument (~ 3 keV). 

Many models of the fast variability in X-ray binaries tie the 
variations in flux to processes in the accretion disc. For example, 
disturbances can be introduc ed at large radii , and be modulated 
as they propagate inwards (iLyubarskiill 19971) . To give rise to a 
QPO, a given mode of variability must modulate the flux. In dif- 
ferent models, the flux modulation is realised in different phys- 
ical regions, and variability arising in the disc may therefore be 
observed at en ergies above those of th e direct disc component. 

Following lAxelsson et al.l (120051) we adopt the relativistic 
precession model. For this model to give rise to QPOs, a given 
precession frequency, and thereby a certain radius, must be pre- 
ferred. In the disc truncation scenario, a natural candidate for 
such a radius is the boundary region between disc and hot in- 
ner flow. Variability transferred from the disc is thus observed as 
modulations in the Comptonized component. The frequency of 
these modulations is determined by the location of the boundary 
region, i.e. the i nner radius of the disc . Adopting the frequency 
identification in lAxelsson et al.l ( 120051) . we can thereby calculate 
a value for the inner radius R m of the accretion disc (in units of 
gravitational radii) using the observed frequencies: 
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where M is the mass of the black hole and v per the perias- 
tron precessional frequency, here identified with the second 
Lorentzian component. This value is for all hard edge observa- 
tions 50 gravitational radii, indicating that the disc is truncated 
far from the innermost stable orbit. We note that this is in agree- 
ment w ith p revious result s by Gilf anov. Churazov & Revnivtsevl 
(1999) and Ibragim ov et al.l (120051) . who find correlations be- 
tween reflection, spectral index and temporal characteristics: as 
the frequencies become lower, the spectrum hardens and reflec- 
tion is reduced. 

The relativistic precession model identifies the two lower fre- 
quencies of the PDS with the nodal and periastron precessional 
frequencies, and a natural next step would be to try and identify 
the third component, or at least tie it to a physical process. 

With the two lower frequencies assumed to arise from rela- 
tivistic precession at the inner edge of the accretion disc, we may 
easily calculate the orbital frequency, v^, as 
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Fig. 4. Correlation between frequencies of the second and third 
Lorentzians. The solid line indicates the orbital frequency (v^) 
for a 1OM black hole as a function of periastron frequency 
(vper). The dashed line is the best-fit power-law to the data, yield- 
ing an index of 0.29 + 0.05, whereas the orbital frequency is pre- 
dicted to follow the relation oc >A£. 



where m is the ma ss of the black hole in units of solar masses 
(IStella et al.| [l999). Figure [4] shows the measured frequency of 
the third Lorentzian component and predicted orbital frequency 
according to Eq.|2]with a black hole mass of 10 M . The dashed 
line in the figure is the best-fit power-law to the data, giving an 
index of 0.29 + 0.05. As can be seen, this is less steep than the 
prediction for the orbital frequency, oc >A*. The frequency of 
the third Lorentzian component thus appears to scale as v° 5 . In 
both the framework of relativistic precession and many other the- 
oretical models for the QPOs, the observed frequencies are ex- 
pected to vary with v^; however, we have found no natural candi- 
date for a frequency proportional to 5 in any model. We stress 
that relativistic precession does no t rule out additional co mpo- 
nents at higher frequencies; indeed. IPottschmidt et al.l(l2003l) re- 
port a frequency component even higher than the ones discussed 
here. 

It is noteworthy that there has been no report of kHz variabil- 
ity in Cyg X-l. If the identification here is correct, this absence 
is naturally explained, as such frequency ranges are far above 
what can be expected from the inner edge of the disc. 

3.2. Hardness-flux correlation 

An interesting feature of the hard edge observations is the be- 
haviour of the hardness-fl ux correlation. As shown by several 
authors (see for exa mple IZhang et ail 119971: IWen et all [l999; 
Zdzi arski et al.l 12002b . this correlation is negative in the hard 
state. This overall trend can be seen in Fig. [2] although there 
is a large spread. 

In order to establish how the hardness-flux relation behaves 
in the hard edge cases, we calculate the Spearman rank correla- 
tion coefficient between hardness (9-20 keV over 2^1 keV) and 
flux for all our individual observations. The result is shown in 
Fig. We find that as Cyg X-l enters the hard edge, the short 
term (~0.5-l hr) negative correlation between hardness and flux 
(within each observation) disappears. We note that the count rate 
remains above 100 counts/s in both channels, indicating that the 
effect is not artificially introduced. 

If we instead combine the observations to study the corre- 
lation on longer timescales, between observations (from several 
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Fig. 5. Spearman rank correlation index between hardness and 
flux as a function of flux. Each point corresponds to one observa- 
tion of ~30 minutes. The hardness is calculated for the 9-20 keV 
over 2—4 keV bands, and the flux is the sum of these two bands. 



hours to days), we find an overall hardness-flux correlation in- 
dex of -0.51 + 0.09. This shows that the correlation on longer 
timescales remains negative, even as it disappears within the 
individual observations. Th e uncertainty is estimated using the 
bootstrap method (see, e.g. JPress et alJl992B . 

As can be seen in Fig. [5] there are a few observations where 
the hardness-flux correlation is actually positive. Looking closer 
at these observations, we find that they were all made within 
a 24-hour period of MJD 50481-50482. Figure |6] shows the 
hardness-flux correlation for one such observation, compared 
with that more commonly seen in the hard state (from an ob- 
servation on MJD 52693). A Spearman rank test of the MJD 
50481 observation yields the value 0.55 + 0.04, showing that 
the correlation is positive. In contrast, the same analysis for the 
normal hard state yields a value of -0.43 ± 0.04. 

To test how our choice of energy bands affects the calcu- 
lated ratio, we extracted lightcurves in the bands matching the 
Standard2 quick-look lightcurves: 2-4 keV, 2-9 keV, 4-9 keV, 
9-20 keV, and 20^-0 keV. This was done for both hard edge 
and typical hard state observations. All possible hardness ratios 
were then calculated. We find that using the 4-9 keV over 2- 
4keV gives a weak correlation for both types of observations, 
probably since the bands are close and narrow. The typical hard 
state observation gave a negative correlation in all cases but one: 
when the highest ranges were used, 20^-0 keV over 9-20 keV, 
the correlation was weakly positive. In the observations from 
MJD 50481-50482, a positive correlation was found for all 
combinations, but varying in strength. 



4. Discussion 

4. 1. The radius of the inner disc 

Spectral studies are often performed on data which do not allow 
precise determination of the black body temperature. In addi- 
tion, spectral modelling of several sources, including Cyg X- 
1, require an additional soft component below 10 keV, mak- 
ing it difficult to determine the true temperature of the inner 
disc. It is therefore difficult to directly determine whether the 
inner radius of the accretion disc differs between spectral states. 
Although the disc truncation model has been successful in ex- 
plaining the full range o f states observed in X-ray binaries 
dDone & G ierlinski 120031) . it is challenged by some observa- 
tional results. Observation of broad iron lines (M iller et al.l2 002: 
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Fig. 6. Correlation between hardness and flux in Cyg X-l in an 
observation from MJD 50481 (upper panel) compared to a typ- 
ical hard state observation from MJD 52 693 (lower panel). The 
correlation in the former observation is weakly positive, in con- 
trast to the typical hard state relation. The hardness is calculated 
for the 9-20 keV over 2-A keV bands, and the flux is the sum of 
these two bands. 



iMiniutti et alll2004t iMiller et al.l|2006l) and a strong black body 
spectral component jMiller et alJT 2006: Rv koff et aLll2007l) both 
seem to indicate an untruncated disc in the hard state. However, 
these interpretations are strongly dependent on the detailed mod- 
elling of the continuum spectral shape, and currently there are no 
observations u nambiguously con flicting with the truncated disc 
model (see e.g. lDone et al.ll2007l) . 

Timing studies offer an independent way of studying vari- 
ations in the accretion flow between different states. Although 
models attempting to describe the PDS are less well established 
than spectral models, many connect observed frequencies with 
Keplerian motion in the accretion disc. Models such as relativis- 
tic precession allow a direct determination of the inner disc ra- 
dius from characteristic frequencies in the PDS. Our results in 
Sect. l3.1l indicate that the disc is truncated far 50 R g ) from the 
innermost stable orbit in the hard edge observations. This sup- 
ports the scenario where the spectrum becomes harder due to the 
inner disc receding, as a result of lower accretion rate. 



4.2. The vanishing hardness-flux correlation 
mini-hysteresis? 



a case of 



Zdzi arski et al.l d2002l) found that the spectral variability of 
Cyg X-l within the hard state on long timescales could be 
described by a combination of two distinct variability pat- 
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terns. In one of the patterns, the spectrum pivots at an energy 
around ~50 keV, creating the observed flux-hardness anticorre- 
lation. In the other, the spectrum varies in normalization with- 
out chang es in the spectral sh ape. The two patterns were mod- 
elled by Zdzi arski et alj (120021) using the Comptonization code 
eqpair dCoppil 119991) . They found that the pivoting pattern 
could be explained by a variable input of soft seed photons to the 
Comptonizing flow. The varying strength of the spectrum with 
constant shape was interpreted as a variable bolometric lumi- 
nosity caused by changes in the local accretion rate and optical 
depth of the hot flow, accompanied by a lower electron temper- 
ature to keep the Compton parameter y = \rkT e /m e .c 1 constant 
(where r is the optical depth and T e the electron temperature 
of the hot flow, and m e the electron mass). Mal zac et alj (2006) 
found that the same two patterns could describe the variability 
on shorter timescales, from hours for the normalization pattern 
to days for the pivoting pattern, in an intermediate state of the 
source. 

Our data show that the anti-correlation between hardness and 
flux resulting from the pivoting behaviour still continues on long 
timescales into the highest hardness levels of Cyg X-l, and sug- 
gest that our hard edge corresponds to an extremely low soft 
photon influx and a maximally truncated disc. The hardness-flux 
correlation within each observation shows that the pivoting be- 
haviour is observable on timescales as short as 0.5-1 hour in the 
normal hard state. The vanishing hardness-flux correlation in the 
hard edge suggests that at this extreme, this variability pattern is 
not observed on such short timescales, even if it is still present 
on long timescales, between different observations. 

A possible interpretation of this behaviour is that when the 
disc is far away, ^50 R g , there is not enough time for the in- 
ner disc radius to adjust to a changing accretion rate within 
the time for a single RXTE observation. A change in the con- 
figuration of the accretion flow and the inner radius of the ac- 
cretion disc is believed to t ake place on the viscous timescale 
dFrank. King & Raindl2003l) . 

(H \ ^ 
-^■J ?dyn , (3) 

where a is the viscosity parameter, H/R the ratio between disc 
scaleheight and radius, and fd yn the dynamical timescale. Our 
frequency analysis suggests that the short term variability both 
in the normal hard state and in the hard edge is consistent with 
a variability of the inner radius of the order of a few R g . If we 
assume that the varying influx of soft photons driving the piv- 
oting behaviour is caused by the variability of the inner radius 
of the accretion disc we can calculate the expected timescale 
for this variability pattern. Assuming the dynamical timescale 
to be the orbital times cale in Fig. [4] and inse rting typical values 
(a = 0.1, H/R = 0.0 1 ; iMaver & Pringldl2007l) . we can estimate 
the timescale of these variations of the inner disc radius to be 
of the order of 0.5-1 hour when the disc is truncated at ~30 R g , 
consistent with the observed timescale for the pivoting pattern in 
the (normal) hard state. 

From Eq. <[3j it is evident that the viscous timescale scales 
directly with the dynamical timescale. Identifying the dynami- 
cal timescale with the orb ital one gives: f v j s c oc t$ oc R^l 2 (for a 
thin disc with constant a: iFrank et al.l l2003). A change in inner 
radius from ~30 to ^50 R g in the hard edge would then result in 
a more than doubled timescale for the same fluctuations in the 
hard edge, compared to the normal hard state. This could then 
explain why we observe the hardness-flux correlation associated 
with the pivoting variability pattern within the time frame of a 



single RXTE/PCA observation in the normal hard state but not 
in the hard edge. Note that the existence of the anticorrelation 
between flux and hardness between different observations, sep- 
arated by several hours to days, in the hard edge proves that the 
same mechanism is present - an increase/decrease of the accre- 
tion rate is accompanied by a reconfiguration of the accretion 
disc with a change of the inner radius. However, the response at 
these large radii is too slow to be observed within the observation 
window of an RXTE/PCA observation. As a result, we may ob- 
serve a change in flux within an observation, but the correspond- 
ing change in hardness only between different observations. 

The effect could be seen as a kind of "mini-hysteresis" and 
compared to that observed in state transitions of certain X-ray 
transients. In e.g. GX 339-4, especially the har d-to-soft transi- 
tion i s observed at a wide range of luminosities (Zdziarski et al. 
2004), seemingly in contradiction to a direct relationship be- 
tween spectral state, luminosity and accretion rate. In some cases 
a luminosity as high as 30 per cent of the Eddington luminosity 
has been observed in the hard state before a transition takes place 
to the soft. It seems as if the local accretion rate increases much 
faster than the accretion flow is able to adapt and change struc- 
ture from a hot flow to a cool disc. Since the source is a transient, 
at the beginning of an outburst the vicinity of the black hole is 
virtually devoid of any material. In this situation the disc requires 
a longer time to build up, and the observed state is thus a function 
of both accretion rate and previous behaviour. In Cyg X- 1 , which 
is a persistent system where the accretion disc is always present, 
the transition between the hard and soft states takes place at the 
same luminosity and no such hysteresis effects are present. The 
effect seen in the hard edge, where we can observe changes in lu- 
minosity that are followed by a delayed change in the inner disc 
radius, can however be seen as a case of "mini-hysteresis" in the 
sense of a delayed response of the disc to a change in accretion 
rate. 

Some of the observations presented here (MJD 50481- 
50482) show a positive hardness-flux correlation within one sin- 
gle observation. This behaviour is the exact opp osite of the nor- 
mal pi voting behaviour. Within the framework of Zdziarski et al.l 
(t2002h . such a variability pattern would be expected if the local 
accretion rate (and optical depth of the hot flow) increased while 
the electron temperature remained unchanged. This could be the 
result of a momentary increase in accretion rate, on timescales 
shorter than that required for a general response in the accretion 
flow. 

5. Conclusions 

We have studied observations of Cyg X-l made during hard-state 
episodes where the hardness ratio is at its highest. In particular, 
the PDS and the hardness-flux correlation have been compared 
to those of the hard state at usual hardness levels. We find sys- 
tematic differences in the hard edge cases: the PDS components 
are found at the lowest end of their frequency range and the usu- 
ally negative hardness-flux correlation weakens or even becomes 
positive within a single observation; however, the overall trend 
between different observations still follows the normal negative 
correlation, consistent with these harder energy spectra occur- 
ring predominantly during low flux episodes. 

We interpret these changes as indications of the inner 
disc being truncated at a large radius, <: 50 R m . The third 
Lorentzian represents a frequency higher than the two generally 
observed, iden t ified a s the periastron and nodal frequencies in 
lAxelsson et al.l d2005l) . This third frequency scales as v" 5 , and 
has no immediate physical interpretation. Our prediction for the 
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orbital frequency also explains why no kHz QPOs are expected List Of Objects 

from Cyg X-l. We relate the timescale for the observed short r-i 

term pivoting behaviour with the viscous timescale for varia- ygnus - on page u 

tions of the inner radius of the accretion disc of a few R g , and 

estimate that this may be more than twice as long in the hard 

edge than in the normal hard state, thus explaining the vanishing 

hardness-flux correlation in the hard edge. Our results support a 

framework where the hard state is explained by a truncated inner 

disc. 
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